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ABSTRACT 

It has been known for a long time that the clustering of galaxies changes as a function of 
galaxy type. This galaxy bias acts as a hindrance to the extraction of cosmological information 
from the galaxy power spectrum or correlation function. Theoretical arguments show that 
a change in the amplitude of the clustering between galaxies and mass on large-scales is 
unavoidable, but cosmological information can be easily extracted from the shape of the power 
spectrum or correlation function if this bias is independent of scale. Scale-dependent bias is 
generally small on large scales, k < 0.1 h Mpc~', but on smaller scales can affect the recovery 
of D. m h from the measured shape of the clustering signal, and have a small effect on the Baryon 
Acoustic Oscillations. In this paper we investigate the transition from scale-independent to 
scale-dependent galaxy bias as a function of galaxy population. We use the Sloan Digital Sky 
Survey DR5 sample to fit various models, which attempt to parametrise the turn-off from 
scale-independent behaviour. For blue galaxies, we find that the strength of the turn-off is 
strongly dependent on galaxy luminosity, with stronger scale-dependent bias on larger scales 
for more luminous galaxies. For red galaxies, the scale-dependence is a weaker function of 
luminosity. Such trends need to be modelled in order to optimally extract the information 
available in future surveys, and can help with the design of such surveys. 
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1 INTRODUCTION 

Galaxies are not expected to form a Poisson sampling of the dis- 
tribution of matter in the Universe. Indeed, it has been known for 



ferent clustering strengths ( Davis & Geller 


19761; iDressleJ 


198C 


ParketalJ 1994; Peacock & Dodds 1994 


Seaborne ct al. 


1999 


Norberg et alj2001, 2002; Zehavi et alj2002 


,2005l;lLi etalj2006) 



showing that they cannot all have a simple relationship linking their 
distribution with that of the matter. This galaxy bias severely lim- 
its our ability to extract cosmological dat a from galaxy surveys 
jPercival et alj|200llsan"chez & Coldl2008t> . 

The large-scale shape of the linear matter power spectrum is 
dependent on Q.„,h because of the chang e in the evolution of the 
Jeans scale after matter-radiation e quality (Silk 1968; Peebles & Yu 
197d;ISunyaev & Zel'dovichll970l ; lBond & Efstathioul 19841 1987k 



Holtzmanlll989ir However, changes in the general shape of the 
power spectrum, such as that caused by this physical process, are 
hard to separate from galaxy bias, which also imprints a signa- 
ture that changes smoothly with scale. Galaxy bias can also affect 
mode-coupling in the power spectrum by changing the non-linear 
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scale, which can le ad to small changes in the Baryon Acoustic Os - 
cillation positions dCrocce & Scoccima rro 2008; Matsubara 20081) . 
Such effects are beyond the scope of our work, and we concentrate 
here on the broad changes to the shape of the galaxy power spec- 
trum. 

The simplest model of galaxy bias is local, linear, determin- 
istic bias, S g (x) = b]j n 6]i n (x) where S g is the galaxy overdensity 
field, and 6] m is the linear mass overdensity field. In this model, the 
bias b lin is constant in space, but can change for different galaxy 
populations. A more generic model may also include a stochas- 
tic element which enters the description as an additional term, 
5 g (x) - £>i in <5i in (x) + 6. Or the bias could be non-local, where it 
depends on a smoothed version of the density field. In this paper, 
we will only be concerned with the relation between the galaxy 
and mass power spectra, and therefore define a practical measure 
of galaxy bias 



PJk) = b(k?PUk). 



(1) 



Such a model can incorporate some of the complexities discussed 
above, but is local in £-space, so cannot include any mode coupling 
terms, which are expected to be present. 

As we move to large scales, the galaxy bias is expected to 
tend towards a constant value . The simplest model for this is the 
peak-background split model jBardeen et ai]|l986l : ICole & Raised 
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1989). Here galaxy formation depends on the local density field. 
Large scale density modes can alter the local galaxy number den- 
sity by pushing pieces of the density field above a critical thresh- 
old. On large scales we expect a linear relationship between the 
large scale mode amplitudes and the change in number density, 
so the shape of the galaxy and mass power spectra are the same. 
Interestingly, such a linear relations hip is broken if the density 
field h as a non-Gaussian component I Dalai et aljiooj j Slosar et"al1 
2008). On small scales we expect galaxy clustering to be different 
from that of the mass, as pai rs of galaxies inside single collapsed 
structures become important ( Seliak 200rj|; |Peacock & Smithl2 000; 
ICoorav & Shethll2002h . 

Although large-scale bias is seen as a nuisance when ex- 
tracting cosmological information from galaxy power spectra, it 
does provide co nstraints on possible galaxy formation models. 
Wild et al. ( 2005) proposed bivariate lognormal models of relative 
bias motivated by observations of galaxy distributions. The data 
were found to support a small (everywhere <5%), but significant, 
amount of stochasticity and non-linearity in these models on all 
scales, implying that galaxy formation is not solely a function of 
local density; these effects must be understood in order to prop- 
erly utilise the next generation of galaxy surveys. A related work, 
IConwav et al.l d2005t) . also supported these findings. 

In order to use the power spectrum to extract cosmological 
information, we need to investigate the transition between scale- 
independent and scale-dependent galaxy bias. This transition is 
known to be a function of th e ga laxy population chos en. Compar- 
ing work in lCole et alj ( 2005 ) and Percival et alJd2007|). which use d 
the same techniques, and the analysis of I Sanchez & Cold \200A 
shows that there are deviations between the shapes of the 2dFGRS 
and SDSS galaxy power spectra. If this results from galaxy bias and 
the fact that the 2dFGRS selected galaxies in a blue band whilst the 
SDSS selected galaxies in a red band, then we should expect that 
similar changes show up in re d and blue galaxies d rawn from a sin- 
gle survey. Using the SDSS, IPercival et alj f2007h showed that at 
k = 0.2 h Mpc~' , the shape of the power spectrum is a strong func- 
tion of the r-band luminosity. In this paper we extend this work by 
splitting in galaxy colour and luminosity, and by fitting models to 
the resulting power spectra, in order to see if the SDSS contains suf- 
ficient galaxies to fully explain the trend observed between SDSS 
and 2dFGRS galaxies. 

Previous work examining the observed depend ence of the 
large scale bias on luminosity was carried out by Norber g" et alj 
d200lh who proposed the p henomenological model b = b\ + 
boLIU, lT7gmark et alj ( f2004l) extended this to better fit the avail- 
able data by including a n extra parameter in t he form b = by + 
btL/Li, + b 3 (M - M*). ISwanson et al. I J2008h reexamined both 
these models for samples split by galaxy colour. IWild et alj (2005) 
also ex amined the colour dependence of relative bias. lZehavi et alj 
J2005I) presented colour dependent and luminosity dependent mea- 
surements of the projected correlation function of SDSS galaxies. 
In all cases subsets of the data existed where some aspect of the 
bias model in question was shown to be a function of colour or 
luminosity. 

We first describe the SDSS catalogues we use in Section [2] 
and how we model the radial selection function of our subsamples 
using luminosity function fits in SectionfS] Section|4]describes our 
calculation of power spectra and uncertainties, and the bias models 
to be fitted are presented in Section [5] Results of fitting the bias 
models are given in Section[6] In Section|7]we present simple mod- 
els for the luminosity dependence of our bias parameters and the 



bin 



absolute 
magnitude 
range 



mean 
Mo.ir 



galaxy 
count 



red 1 
red 2 
red 3 
red 4 
red 5 
red 6 

blue 1 
blue 2 
blue 3 
blue 4 
blue 5 
blue 6 



-22.30 < M 0Ar < -21.35 



-21.35 < M 0Ar < 
-20.89 < M 0Ar < 
-20.47 < M 0Ar < 
-20.00 < M 0Ar < 
-19.34 < M 0Ar < 

-22.30 < Mo.!,. < 
-21.35 < M 0Ar < 
-20.89 < Mo. 1; . < 
-20.47 < M 0Ar < 
-20.00 < M 0Al . < 
-19.34 < M 0Ar < 



■20.89 
•20.47 
■20.00 
■19.34 
•17.00 

■21.35 
■20.89 
■20.47 
■20.00 
■19.34 
■17.00 



-21.79 ±0.27 
-21.12 ± 0.13 
-20.70 ±0.12 
-20.27 ±0.14 
-19.75 ±0.19 
-19.01 ±0.38 

-21.67 ±0.24 
-21.22 ±0.13 
-20.69 ±0.12 
-20.27 ±0.14 
-19.74 ±0.19 
-18.91 ±0.49 



49167 
41462 
37819 
34651 
29742 
17582 

17480 
25208 
28928 
32066 
36889 
48754 



Table 1. Description of the subcatalogues analysed in this paper. The lim- 
its of the absolute magnitude bins are given, with the weighted mean and 
standard deviation. We also give the galaxy count for each bin. 



results of fitting these to our data. Conclusions and discussion are 
presented in Section[8] 



2 GALAXY CATALOGUES 

We use galaxy catalogues selected from the Sloan Digital Sky 
Survey DR5 main galaxy sample. The Sloan Digital Sky Survey 
(SDSS: lYork et aljboOOr , lAdelman-McCarthv et alj|2006l), which 
was recently completed, used a 2.5m telescope dGunn et alJl2006f) 
to obtain 10 4 square degrees of imaging data in fi ve passbands u, g, 
r, i and z (Fukugita et al |l99dlGunnetalJl998h . The main galaxy 
sample JStrauss et alJl2002h consists of galaxies with Petrosian r- 
band magnitude m rpe , rosian < 17.77. This gives approximately 90 
galaxies per square degree, with a median redshift z = 0. 1 1 ; in thi s 
paper we use the DR5 sample ( Adelman-McCarthy et al. 2006). 
We exclude a small subset of the data taken during initial sur- 
vey operation, for which the apparent magnitude limit fluctuated. 



This gives 410 095 galaxies with 14.5 < m r . 



< 17.77 and 



redshift > 0.003; the lower redshift cut strongly reduces the con- 
tribution from mis-c lassified stars. This is the catalogue used in 
IPercival et al. and further details can be found here. 

Where specified, we have K- corrected the galaxy l uminosities 
using the methodology outlined in iBlanton et ail < l2003 jlbh . We also 
use the same z = 0.1 shifted r-band fil ter to define our luminosities 
(as discussed in Blanton et al. 2003b), which we refer to as M0.1,. 
throughout this paper. Absolute magnitudes and fc-corrections were 
calculated assuming h = lOOkms-'Mpc" 1 , fl M = 0.3 and fi A = 
0.7, and we have applied the reco mmended AB corr ections to the 
observed SDSS magnitude system dSmith et alj|2002l) . 

This sample of main galaxies was split into 6 absolute magni- 
tude limited subsamples giving approximately equal numbers of 
galaxies in each bin. Each subsample was then further divided 
into red and blue subsamples defined by a constant colour cut of 
Mq.i ~~ M Q lr = 0.8. Fig. [T] illustrates our chosen simple colour 
split, showing that this cut clearly divides the two populations. We 
have tried more complicated cuts without significant change in our 
results. 

The absolute magnitude cuts and number of galaxies in each 
of these subcatalogues is given in Table[T] Fig.f2]illustrates the dis- 
tribution of red and blue galaxies within the imposed absolute mag- 
nitude limits. At a given absolute magnitude the high and low red- 
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M 01r -5log 10 (h) 

Figure 1. The distribution of galaxies in colour-absolute magnitude space. 
The dashed line denotes the boundary of the colour split at Afn.l s - M).l r = 
0.8. 

shift limits result from the apparent magnitude limits of the survey, 
this can be seen in Fig. [2] We could have cut the catalogues in order 
to make these sub-catalogues volume limited, but we wish to retain 
as much signal as possible. Significant cuts would have been re- 
quired in order to rem ove all effects of evolution and K-corrections. 
|Percivaletal]j2007h showed a compromise using "pseudo-volume 
limited" catalogues. We choose instead to estimate the redshift dis- 
tribution by fitting the luminosity function (see the next section), so 
we can model an apparent magnitude cut as easily as a cut in abso- 
lute magnitude. This approach allows us to work with considerably 
more galaxies and retain the maximum information. 



3 MODELLING THE RADIAL SELECTION FUNCTIONS 

For each of the catalogues described in the previous section, we 
need to model the radial selection function. The angular mask is 
the same for all catalogues as all the cuts being appli ed are inde- 
pende nt of angular position. We use a HEALPIX0(G6rskietal. 
2005) mask to describe the angular galaxy distribution as described 
in lPercival et al. We derive the radial distribution of galax- 

ies in each of our samples from fits to the luminosity functions for 
either red or blue galaxies as appropriate. Section [3Tl describes our 
luminosity function model and Section [T2l describes our method 
for transforming this into a redshift distribution. 

3.1 Redshift evolution corrected luminosity functions 

Our subsamples of red and blue galaxies contain sufficiently large 
numbers of galaxies to allow us to calculate redshift evolution 
corrected luminosit y functions. We have found that well known 
Schechter function dSchechte n [l976l) fits the data well. In terms 

1 http://healpix.jpl.nasa.gov 
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Figure 2. The distribution of galaxies in absolute magnitude-redshift space, 
red galaxies are shown in the top panel, blue galaxies in the lower panel. The 
horizontal, dashed lines denote the boundaries of the absolute magnitude 
bins. 

of absolute magnitude and with modifications to include redshift 
evolution the Schechter function is given by 

$>(M,z)dMdV = »0.41og e (10) 10 m ^o)p 

l0 0A(M*-M-Qi : - ZB m l+a)exp{ _ lQ 0A(M*-M-Q( z - Z0))]dMd ^ (2) 

where P an d Q are redshift evolution parameters following the con- 
vention of iLin et alj J19991) : P allows for density evolution, while 
Q allows for luminosity evolution. We have also tried fitting non- 
parametric models, models with so many parameters that the final 
shape of the function is to an extent indepe ndent of the shapes o f 
the contributing terms, based on the work o f lBlanton et aIU2003bl) . 
However we find no significant change in the resulting redshift dis- 
tribution, and we therefore only consider fitting a Schechter func- 
tion in the remainder of our paper. 

The luminosity function parameters were determined using a 
maximum likelihood method dLinetal.l l 1999) implemented with 
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Figure 3. The best fit redshift-evolution corrected Schechter functions 
(Eq[2] see text for details) representing the luminosity functions of the red 
galaxies, shown as the solid red line, and blue galaxies, shown as the dashed 
blue line, both at a redshift of z = 0.1. The red, vertical crosses and blue, 
diagonal crosses represent the 1/V„ MJ data estimates of the true luminosity 
function for red and blue galaxies respectively, redshift evolution corrected 
using the respective best fit parameters for Eq.ffjto a redshift of z = 0. 1. 



the minimisation routine Powell (Pre ss et al .Il992h . Convergence to 
the Likelihood maximum was confirmed by starting the minimisa- 
tion routine at a number of widely separated initial parameter sets, 
and observing that the same best-fit parameters were obtained. The 
results for the red and blue galaxy samples can be seen in Fig. [3] 
Note that the overall normalisation, h cannot be determined by this 
maximum likelihood method and we determine it directly from the 
data. Further discussion of our procedure to calculate luminosity 
fu nctions and the parame ters of the resulting fits will be presented 
in lCresswell et al.l d2008h . In this paper we simply use the luminos- 
ity functions to estimate redshift distributions for the selected cat- 
alogues, so the recovered parameters are not important, provided 
the fits are an adequate match to the data: for instance, there is 
some apparent discrepancy between the l/V mov data estimates and 
the best fit Schechter function for the red galaxies near an absolute 
magnitude of M^i r w -19.3; as the luminosity function is slowly 
changing in this region and the normalisation is determined inde- 
pendently of the fit this will have no significant effect on the shape 
of any derived redshift distributions, which would here be domi- 
nated by the effect of the apparent magnitude limits of the survey. 



3.2 Redshift distributions 

Given a luminosity function as defined in Eq. l[2} we can integrate 
to determine the redshift distribution using 



J 



upper rfy 

®(M',z)dM' — dz, 



(3) 



where Mi owcr is the minimum of the lower survey absolute magni- 
tude limit at a given redshift, and the lower bin magnitude limit, 
and M U p pcl is similarly defined from the upper boundaries of the 
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Figure 4. Histogram showing the redshift distribution of red galaxies in 
absolute magnitude bin 2. The thick black line is the model distribution 
derived from the red galaxy luminosity function via the relation shown in 
Eq. 0). 



sample and bin. An example of the data and model redshift dis- 
tributions for one of the subsamples, the red galaxies in absolute 
magnitude bin 3, is shown in Fig. [4] Good agreement is seen in this 
plot, and between model and data redshift distribution for all of our 
subsamples, validating our procedure. 



CALCULATING POWER SPECTRA AND 
UNCERTAINTIES 



Power spectra were calculated as descri bed inlPercival et al. 

Ij2007h . 

using the standard Fourier technique of iFeldman et alj d 19941) . The 
radial and angular selection functions were included in this method 
by creating a random catalogue with the same selection func- 
tion as the galaxies, but with Poisson sampling and lOx as many 
ga laxies. The galaxies were weighted using the optimal weights 
of lFeldmanetal.Ul994h . Obviou sly, no additional bi as-dependent 
weighting scheme (such as that in lPercival et al.l20 04 was applied). 

For the SDSS DR5 sample we have cre ated 2000 Log-No rmal 
(LN) catalogues, calculated as described in ICole et alj J2005h . We 
assumed a flat ACDM power spectrum with Sl„,h = 0.2, and 
£l;,/fl„i = 0.15. Normalisation was matched to that of L, galaxies 
as defined below. We have applied a colour and luminosity depen- 
dent bias model that is scale-independent to these catalogues. The 
scheme was iteratively mat ched to our results: in itially we used the 
luminosity-bias relation of iNorberg et alj d200ll) to calculate cata- 
logues and the corresponding covariance matrix for the data. Hav- 
ing fitted the data, we then calculated new catalogues with a colour- 
luminosity bias scheme that was a better fit to the data. This change 
resulted in a negligible change to the best-fit bias model, so we are 
confident that our results do not depend on this choice. 

Power spectra were calculated from the mock catalogues using 
exactly the same process as for the actual data, and the covariance 
matrices were estimated for each of the subcatalogues described 
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in Section [2] Because LN catalogues for each subcatalogue were 
drawn from the same underlying density fields, we use these mocks 
to calculate correlations between the power spectra for different 
subcatalogues: these will be correlated as the volumes overlap. 



5 FITTING MODELS OF GALAXY BIAS 

As discussed in Section [T] we define bias as the ratio between 
the galaxy power spectra and the linear matter power spectrum 
(Eq.[T). We fit the observed power spectra with models calculated 
from a linear model p ower spectra, using the fitting formulae of 
lEisenstein & HvJ ( |l998h . with parameters given by the concordance 
cosmology Cl = 0.241, fl A = 0.759, H = 73 .2, cr g = 0.761, 
n s = 0.958, Q. b /Sl m = 0.175 dSpergel et ai]|2006l) . multiplied by a 
bias model. The model is then conv olved with the window function 
for each catalogue, as described in IPercival et af] J2007h . We will 
be interested in relatively large scales k < 0.4/iMpc~' and assume 
that, on these scales, the power spectrum band-powers result from 
a multi-variate Gaussian distribution, and that they are correlated 
within a single power spectrum, and between different power spec- 
tra. We perform Maximum Likelihood fits to the measured power 
spectra for three models of the bias b(k). First we assume that the 
bias is constant and fit on scales k < 0.21 /iMpc~'. We then con- 
sider two scale-dependent models for b(k), fitting to smaller scales. 

Results are given in section[6] 

First we use the g-model JCole et alj2005T) . 



b(k) 



1 + Qk 2 
l+Ak' 



(4) 



is the asymptoti c large-scale bias, and Q and A are pa- 



Cole etalJ ( [2005T) used halo model catalogues to suggest 



where by,, 
rameters. 

that A = 1.4 in redshift-space leaving just a single nuisance pa- 
rameter Q to be fitted to data. Q would be expected to change for 
different galaxy populations, and it is this that we test as a function 
of galaxy colour and luminosity. 
We also consider a model with 



b(k) = h 



1 + ■ 



bf Pu 

lin " 



(5) 



where the parameter P acts as an additional shot noise term. This 
model has a physical basis as this term could account for a change 
in the shot noise. This would arise if halos Poisson sampling the 
density field, and galaxies are located in those halos, i.e. P is the 
contribution t o the galaxy power spec t rum from the one-h alo term 
ISeliakl i l200lh : ISchulz & Whit3 j2006h : lGuzik et aljj2007h . We will 
refer to this model as the P-model. 

The relative merits of these two models are discussed in 
who argue in favour of the ease of physical in- 
terpretation of the P-model. In addition, although they show that for 
classes of cosmological model containing a free-streaming hot dark 
matter component consisting of relic thermal axions, the g-model 
becomes highly pathological due to a degeneracy between Q and 
the particle mass. They suggest that this pathology may extend to 
other models with light thermal relic components. 



6 RESULTS FROM FITTING THE POWER SPECTRA 

We have fitted bias models to power spectra calculated from the 
catalogues parametrised in TableQ] First we consider fitting scale- 
independent bias on large scales, and then we consider scale- 
dependent bias. All bias measurements were calculated relative to 
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Figure 5. Comparison of relative, large-scale, constant, linear bias as a 
function of luminosity for galaxies split into red and blue colours. Five 
and six pointed stars represent our measurements of bias relative to our red 
galaxy M+ bin, for r ed and blue galaxies re spectively. Squares and circles 
show the results from Swanson et af |2008) renormalised to match our M* 
values, see text for details. The upper, solid line is a fit to Eq.|6]for our red 
galaxy bias points, see text for details. The lower, dash-dot line is the fit to 
blue galaxies. The horizontal dashed line shows no bias relative to our red 
M„ bin. 



the bias of red galaxies of luminosity L*, fc* Jc d, where was cal- 
culated from the fit to the red-galaxy luminosity function. 



6.1 Large scale bias 

If we assume that the bias does not change with scale, and fit to very 
large scales k < 0.21 hMpc~\ then the resulting bias amplitudes, 
measured relative to fc», re d, are shown in Fi g. [5] For comparison we 
also plot the data of ISwanson et alj d2008l) with an average evolu- 
tion correction removed from each magnitu de bin. In order to com- 
pare samples with different values the ISwanson et al ] J2008h 
data have been offset to so that the linear interpolation of the data 
points either side of our L* bin passes through our bin, these cor- 
rections are small, as -0.05 Mo.i, for red galaxies and as +0.01 Mo.i,. 
for blue galaxies. As can be seen, we recover the same trends with 
the bias of blue galaxies monotonically increasing with galaxy lu- 
minosity, while the red galaxies show more complicated behaviour 
with increased bias for both high and low luminosity galaxies. This 
comparison is discussed further in Section[8] 



6.2 Scale dependent bias 

The best-fit parameters for the Q-model and P-model fits are pre- 
sented in Table|2] These numbers were calculated by fitting each of 
the 12 power spectra with the bias model multiplied by the linear 
power spectrum. The best-fit values and their trends as a function 
of galaxy colour and luminosity are analysed further in Section|7] 

The model and data power spectra are compared in Figs. [6] 
&|7] These data are divided by our fiducial model convolved with 
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Figure 6. Comparison of best-fit power spectra calculated with bias from the Q-model (solid lines) to the data (circles with lcr errors). Each power spectrum 
is divided by our fiducial linear power spectrum convolved with the appropriate window function. 
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Figure 7. As Fig. [6] but now for model power spectra calculated assuming 
the P-model for galaxy bias. 



Q-model P-model 



bin 




/ 


Hin 


Q 




bun 


P 


red 1 


1 


40 


±0.02 


9.45 ± 0.70 


1 


35 ± 0.02 


512 ±48 


red 2 


1 


22 


±0.03 


9.24 ± 0.83 


1 


17 ± 0.03 


375 ± 40 


red 3 


1 


17 


±0.03 


9.44 ± 0.94 


1 


13 ± 0.03 


351 ±42 


red 4 


1 


20 


±0.04 


7.72 ±0.93 


1 


16 ± 0.04 


277 ± 46 


red 5 


1 


25 


±0.05 


7.26 ± 1.03 


1 


20 ± 0.05 


272 ±51 


red 6 


1 


44 


±0.08 


8.03 ± 1.45 


1 


39 ± 0.07 


419 ±80 


blue 1 


1 


09 


±0.03 


13.74 ± 1.29 


1 


04 ± 0.03 


541 ±48 


blue 2 





96 


±0.03 


9.89 ± 1.34 





92 ± 0.03 


274 ± 40 


blue 3 





94 


±0.04 


7.74 ± 1.43 





90 ± 0.04 


177 ±43 


blue 4 





89 


±0.05 


7.44 ± 1.74 





86 ± 0.05 


151 ±46 


blue 5 





91 


±0.07 


4.64 ± 1.78 





87 ± 0.06 


66 ±51 


blue 6 





92 


±0.14 


2.99 ± 2.97 





87 ±0.13 


17 ±80 



Table 2. Maximum Likelihood parameters calculated by fitting the Q-model 
and P-model for scale-dependent galaxy bias as given in Eqns.|4]&[5] to 
power spectra calculated for the catalogues described in Section [2] Best-fit 
values are presented together with 1-<t uncertainties. 

the appropriate window function. Comparisons are shown for the 
red galaxies in the left hand panels and blue galaxies on the right, 
with the rows of panels showing different absolute magnitude bins 
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with the brightest at the top and faintest at the bottom. The offset 
visible between the data and model in the lower left panel is due 
to the high best-fit large-scale bias in the best-fit solution (Eq.(6] fit 
parameters given in TableO. Such effects arise because the data are 
correlated, and the maximum Likelihood solution does not match 
the ^-by-eye expectation. 



7 MODELLING BIAS AS A FUNCTION OF 
LUMINOSITY 

In this section we try to fit a simple model for the luminosity de- 
pendence of the parameters in the bias models given in Eqns.[4]&[5] 
The asymptotic l arge-scale bias fc lin is known to be a functio n 
of galaxy luminos ity jNorberg et al J 1200 ll ; iTegmark et al.ll2004) . 
and galaxy colour dSwanson et al.l2008h. Here, w e extend the form 
of the model introduced by Nor berg et al.l 1 120010 in order to cope 
with the more complicated behaviour seen when we also split 
galaxies by colour. We assume the three parameter model 



L L+ 
b Ym (L) = flj + a 2 — + a 3 — - 



(6) 



For relative biases where b lm (L+) = 1, we would have that a l + a 2 + 
ci3 = 1, so the model would only have two free parameters. This 
model reduces to the form of iNorberg et al.l j200lh in the case of 
a 3 = 0. 

In addition, we parametrise 



and 



Q(L) = qi + q 2 (M - M*), 



P(L) = p l+ p 2 (M-M*), 



(7) 



(8) 



where q t , q 2 , p t and p 2 are parameters that we can fit to the data. 
We therefore have 5 free parameters in our fit to the data in addi- 
tion to choice of whether to use the P or g-model. We have fitted 
these parameters by performing a simultaneous Maximum Likeli- 
hood search using all 12 power spectra, and allowing for covari- 
ances between band powers in each, and between different power 
spectra. I.e. we do not simply fit the recovered data values in Fig. [8] 
but instead perform a new fit to the data. The resulting bias models 
are compared with the results calculated when we allowed by m , Q 
and P to match each catalogue individually in Fig. [8] As can be 
seen, this simple model does very well in matching the luminosity- 
dependent trends observed. 

The resulting best-fit parameters allow us to define the fol- 
lowing models for redshift-space galaxy power spectra. Using the 
g-model, for red galaxies, 



/\ vfi „<(/-./,-> = jo.92 + 0.10^ +0.19^- 



xP lm (L,k) 
and for blue galaxies, 



1 + (7.5 - 1.0(M - M*))k 2 
1 + Ak 



(9) 



P,. llM ,jL.k) = (0.81+0.07-^+0.02^ 



XPy m (L,k) 



\ + {63-~iA{M-M±))k 1 



1 +Ak 



(10) 



These model power spectra are compared with those observed in 
Fig. [9] where good agreement is seen. 
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Figure 9. As Fig. [6] but now showing models with parameters calculated 
from a simple fit to the luminosity dependent trends observed in Fig. [8] 
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Figure 10. As Fig. |7] but now showing models with parameters calculated 
from our simple fit to the luminosity dependent trends observed in Fig. [8] 



Using the P-model, for red galaxies, 

Pgalred(.L, k) 



0.89 + 0.10— +0.18^- 



xP lto (L,fe) + (200 - 40(M-MJ), (11) 



and for blue galaxies, 



ga{ blue 



(L,k)- 



0.77 + 0.07— +0.02 — 
L* L 



x P hn (L, k) + (160 - 130(M - M*)), (12) 



where, as previously noted in the text, A = 1.4. These model power 
spectra are compared with those observed in Fig. [To] where, as in 
Fig- 13 good agreement is seen. 

The x 1 values of the best-fit models including the P-model 
of bias are good, with 609 given 571 degrees-of-freedom for the 
red galaxies and 439 for the blue galaxies given the same number 
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Figure 8. Top row: best-fit £>/,„, and Q model parameters as a function of absolute magnitude, as presented in Table|2] from a fit including the Q-model bias 
prescription. Bottom row: best-fit &/,„, and P model parameters as a function of absolute magnitude, as presented in Table|5] from a fit including the P-model 
bias prescription. Circles and crosses are for red and blue galaxies respectively. Solid, horizontal lines about each data point show the extent of each absolute 
magnitude bin. Solid lines show the models of Eqns.l9l& ll H and dashed lines show the models of I10l &ll2lin the appropriate panels for red and blue galaxies 
respectively. 



of degrees of freedom for the blue galaxies. For the Q-model, the 
corresponding numbers are 624, and 440. These numbers depend 
strongly on the covariance model adopted and, as a consequence, 
we do not analyse these further other than to note that the fits seem 
to give reasonable numbers. 



8 CONCLUSIONS 

We have quantified the bias of red and blue galaxies as a function of 
luminosity on scales k < 0.4/iMpc~'. We find clear differences in 
the large-scale asympt otic bias between blu e and red galaxies, sim - 
ilar to that found bv lZehavi etal] d2005h and lSwanson et al.l (2008), 
and shown in Fig. [8] At large scales red galaxies have are more bi- 



ased than blue galaxies for all luminosities. The bias of blue galax- 
ies is a strong, monotonically increasing function of luminosity, 
with more luminous galaxies being the most biased. For red galax- 
ies, the picture is mo re complicate d, with an increase in the bias 
of faint red galaxies I Zeha viet al.l J2005h find a similar result but 
stress the scale dependence of the observation). This trend is sig- 
nificant at the 7.2-cr level.We fi nd a stronger bias for lu minous blue 
galaxies, compared with that of lSwanson et al. I j2008l) . The reason 
for this is unknown, althou gh it is worth noting that our catalogue 
is larger than that used bv lSwanson et ID J2008h . as they used ad- 
ditional cuts, constructing volume limited subsamples, compared 
with our sample selection procedure (see section[2}- We performed 
a simple test of the dependence of this difference on observed 
scale, repeating our analysis for scales 0.038 < k < 0.070 h Mpc~' 
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which approximates the lSwanson et alj d2008h measurement scale 
of ~ 20h~ [ Mpc, the discrepancy in the bias of the brightest blue 
galaxies was not removed. However, the change of scales brought 
the best fit red galaxy line ar bias model into better agreement with 
the ISwanson et alj J2008h faint, red galaxies; the model is poorly 
constrained at the faint end due to the lack of unique data points in 
that range. It is worth noting that, due to dif ferent colour cut crite- 
ria, the brightest absolute magnitude bins in lSwanson et al.1 1 2008) 
are inherently more red than ours, this might lead one to expect 
larger bias measurements for those bins, the opposite of the ob- 
served trend. The bright blue sample contains the smallest numbers 
of galaxies, so the errors on the relative bias should be the largest of 
any sample. However, our expected errors are insufficient to fully 
explain the discrepancy, and there remains no obvious reason for 
this difference. 

The work on constant bias models has been extended by con- 
sidering how the turn-off from constant large-scale bias depends on 
galaxy colour and luminosity. We have compared two models for 
this turn-off: the P -model and Q-model given in Eqns. [4]&[5] re- 
spectively. We find that there is little to choose between the two in 
terms of how well they can fit the current data, and both provide 
adequate fits to the power spectrum trends observed as a function 
of galaxy luminosity and colour at the current level of data preci- 
sion. A lthough there is no observational motivation, Hama nn et al.1 
(2008) argued that the P-model has a physical motivation, and 
therefore offers a more attractive model of bias. 

We use the values of P and Q obtained for these two models to 
quantify the degree of divergence from a constant bias model. Note 
that these parameters change both the position and amplitude of the 
turn-off. We find that the best-fit values of P and Q, shown in Fig.[8] 
are a strong function of luminosity for blue galaxies. Red galaxies 
show far weaker evolution with luminosity, and are consistent with 
the hypothesis of no change in the scale at which the bias can no 
longer be described as a constant, to current data precision. Inter- 
estingly, this trend in the turn-off from apparent scale-independent 
behaviour does not match that of the amplitude of the large-scale 
bias. If amplitude and turn-off scale were linked, we might have ex- 
pected the Q and P parameters to match for red galaxies of high or 
low luminosity, where the constant bias component matches, but be 
different for intermediate luminosity galaxies. We do not observe 
such a trend, suggesting that there might not be a simple link be- 
tween the amplitude of bias and the scale at which the 1-halo term 
becomes important. 

It is clear that the simple P and Q-models will become insuf- 
ficient to model the observational data, both on very small scales, 
and as the data improve. They are, after all, simply motivated to 
fit observed trends, and do not encompass all of the physics in- 
volved. One alternative and m ore compl e x pres cription for bias 
prediction is that presented in lYoo et all 1.2008) which can pre- 
dict P bs(k) given a set of cosmological parameters and a HOD 
model, the parameters of which have been determined from the 
observed, projected correlation function (the correlation function 
provides additional information to that of the power spectrum as 
they are sensitive to different scales). However, the potential limits 
of simple, phenomenological models, such as the P and Q-models, 
does not stop us being able to use them to investigate trends in the 
data, as we have done in this paper. The current data clearly show 
that, on relatively large scales (our small scale limit of 0.4 h Mpc~' 
approximately corresponds to 17 /T 1 Mpc), blue and red galaxies 
have very different bias properties. It is the blue galaxies that are 
faint in the r-band whose relative clustering on different scales most 
closely resembles that of the mass, and the luminous blue galaxies 



the least. This supports the hypothesis that the s hapes of the 2dF- 
GR S and SDSS main ga l axy power spectra of ICole et alj fc005h 
and lTegmark et al.1 j2004h : |Percival et al.1 b007h differ because the 
average galaxy bias of each sample differs, caused by different sam- 
ple selections. The clustering of the 2dFGRS galaxies, on average, 
would be expected to be a better tracer of the linear clustering sig- 
nal out to smaller scales than the SDSS galaxies. 

This highlights the importance of sample selection for fu- 
ture galaxy surveys, and the importance of understanding galaxy 
bias for extracting cosmological information from power spectrum 
shapes from such surveys. Red and blue galaxies show very differ- 
ent trends in their large-scale bias and the scale-dependent smaller- 
scale bias as a function of luminosity. It is clear that galaxies need 
to be split into sub-populations by more than just the luminosity in 
a single band in order to properly understand and model bias. 
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